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We have measured the evolution of the tunneling density of states (DOS) in superconduc- 
tor/ferromagnet (S/F) bilayers with increasing F-layer thickness, where F in our experiment is 
the strong ferromagnet Ni. As a function of increasing Ni thickness, we detect multiple oscillations 
in the DOS at the Fermi energy from differential conductance measurements. The features in the 
DOS associated with the proximity effect change from normal to inverted twice as the Ni thickness 
increases from 1 to 5 nm. 

PACS numbers: 74.50.+r, 74.45,+c, 75.70.Cn 



Hybrid systems consisting of superconducting (S) and 
ferromagnetic (F) materials have attracted substantial 
attention due to their interesting properties and poten- 
tial for applications. 1 The superconducting proximity ef- 
fect in such systems is normally short-ranged, due to the 
large exchange energy in the F material. When a Cooper 
pair crosses the S/F interface, the spin- up and spin-down 
electrons enter into different spin bands, and the center 
of mass coordinate picks up an oscillatory factor. 2 The 
physical manifestations of this oscillation can be observed 
as a series of transitions between "0" and 'V states in 
S/F/S Josephson junctions as a function of increasing 
F-layer thickness^ or as oscillations between "normal" 
and "inverted" proximity features in the tunneling den- 
sity of states (DOS) of S/F/I/N tunnel junctions. 5 (Here 
I is an insulator and N is a normal metal.) 

With substantial experimental effort in S/F/S Joseph- 
son junctions, the 0-ir transition has been confirmed by 
many experimental groups £d&r— What is surprising is 
that, unlike in the Josephson geometry, the oscillatory 
behavior of the DOS in S/F/I/N structures has been ob- 
served convincingly only once, as a single normal-inverted 
transition in samples with a weakly-ferromagnetic alloy 
for F£ In experiments using strong ferromagnet s, the re- 
sults have been less clear j 13 i 14 At this time, to the best of 
our knowledge, there is no definitive experimental answer 
to the question of whether the DOS in S/F/I/N struc- 
tures oscillates as a function of F-layer thickness when F 
is a strong ferromagnet. The primary goal of this Rapid 
Communication is to answer this question. 

Our S/F/I/N tunnel junctions are fabricated by ther- 
mal evaporation and sputtering, using a series of me- 
chanical masks (Fig. 1). We first evaporate a 150 
nm strip of Al (N), then we immediately backfill the 
chamber with 300 Torr of a 10% 2 , 90% Ar mix- 
ture. Exposing the freshly evaporated Al to the O2 
quickly (while the Al is still hot) provides good con- 
ditions for oxide growth. The O2 exposure continues 
for w 12 h to produce a robust layer of AI2O3 (I) on 
the Al surface. Next, we change masks and evaporate 
a thick layer (200 nm) of SiO^ to define the junction 
geometry. When using mechanical masks for the top 
leads, shadow effects can cause unwanted regions at the 
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FIG. 1: (Color online). (Left-hand side) Vertical represen- 
tation of the top lead (thin film multilayer) with associated 
thicknesses. The functions of the various layers are described 
in the text. The only variable thickness is Ni, djvi- (Right- 
hand side) Top view of S/F/I/N junction using an optical 
microscope. The junction area A is defined by the edges of 
the SiOx and the AI/AI2O3 bottom lead. The vertical and 
horizontal arrows illustrate the junction width and length. 



edges of the junctions where the Ni thickness is not 
well defined. The SiO^ is in place to avoid the appear- 
ance of edge effects in our data. Finally, we sputter a 
Cu(5nm) /Ni(x) /Cu(10nm) /Nb(150nm) /Au(15nm) mul- 
tilayer. The choice of Ni and Cu is beneficial because our 
Ni has a relatively long spin-diffusion length, 21 ± 2 nm, 
as compared to our maximum Ni thickness, along with a 
low resistivity pNi = 33 ±3 n^m. 15 Ni also provides weak 
asymmetry of spin-dependent scattering in the bulk and 
at Ni/Cu interfaces, and a low average Ni/Cu interface 
resistance. 15 These attributes should simplify theoretical 
analysis. The Cu layer adjacent to the AI2O3 has been 
found to increase the effectiveness of the tunnel barrier. 
The Au deters oxidation of the Nb layer. Throughout 
the process we must break vacuum, but the consistency 
and reproducibility of our results suggest that this has 
little effect on the quality of our junctions. Due to a high 
level of oxygenation of our Al, its T c w 1.9 K. Thus, we 
performed our measurements at 2.1 K, with the Al in the 
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normal state. 

Using a four-terminal lock-in technique, we measure 
the voltage-dependent differential conductance of our 
samples, dI/dV(V), which approximates the DOS of our 
S/F bilayer. (The true DOS would be attained if mea- 
sured at T = K when there is no zero-bias anomaly; 
we will use "DOS" to refer to our non-ideal differential 
conductance measurements.) We normalize dl/dV by 
multiplying by the normal-state resistance Rn, deter- 
mined from the inverse of the differential conductance at 
large bias voltages. Our junction area-resistance prod- 
ucts, ARn ~ 2 x 10 5 Qfim 2 , are substantially higher 
than those of other groups using passive oxidation / 6 ! 17 
but much lower than those produced by specialized oxi- 
dation techniques £^ From the data, it can be seen that 
this has little to no effect on our measurements. 




FIG. 2: (Color online). The blue dot-dashed curve shows 
dl/dV vs. V for a djvi = 4.5 nm junction, where the data 
are slightly smoothed and normalized. The blue dotted line 
is a linear fit to large bias voltage for the previous blue dot- 
dashed curve. The red solid curve shows the same dl/dV data 
after subtracting off the linear background feature. Inset: 
Normalized dl/dV vs. V for an N/I/S junction with a 15-nm 
Cu buffer layer between S and I. All data were taken at 2.1 
K. 



observed in S/F/I/N tunnel junctions measured by other 
groups j 18 i 19 

The inset of Fig. 2 shows a plot of dl/dV for an N/I/S 
junction with a 15-nm Cu buffer layer between the Nb 
and AI2O3 layers. The standard N/I/S junction behavior 
illustrates the quality of the insulating barrier. The Cu 
buffer layer thickness in the N/I/S sample is the sum of 
the Cu layer thicknesses in the S/F/I/N samples. This 
was chosen to illustrate the minor effect of the Cu in 
the S/F/I/N samples as a whole. Even with this rather 
thick layer of Cu and the elevated temperature, we see 
a gap where dl/dV nearly goes to zero. Another sign of 
junction quality is that dl/dV is featureless for \V\ ^> A 
up to 10 mV (not shown), the maximum measured \V\. 
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Our dl/dV data exhibit a small negatively sloped, lin- 
ear background, which appears to be a component of the 
normal state of our junctions. Figure 2 shows a plot of 
dl/dV from our 4.5-nm sample before and after subtrac- 
tion of the linear background, along with the linear fit 
to dl/dV for \V\ > A, where A w 1.4 mV is the gap 
parameter for Nb. There is another normal state charac- 
teristic which we do not correct for in our measurements. 
A slight V-shaped feature centered at V becomes ap- 
parent at large Ni thicknesses. (It is visible on the red 
solid curve of Fig. 2 at \V\ > 2.5 mV.) We did not have 
a magnet on our apparatus to force the top Nb layer into 
the normal state; nevertheless, we emphasize that the 
slight negative slope and V-shaped feature in our back- 
ground are both much smaller than background features 



FIG. 3: (Color online). Normalized differential conductance 
vs. voltage for S/F/I/N junctions with several different Ni 
thicknesses. The Ni thickness, dNi, is labeled under each 
curve at the right. As one moves up the figure, each panel 
has an increasing expansion of the vertical scale where only 
the lowest trace is normalized to 1 and the others are displaced 
upward for clarity. All data are taken at T — 2.1 K. 

Figure 3 shows a plot of the DOS of our tunnel junc- 
tions while varying the Ni thickness. In a separate experi- 
ment at T = 5 K, we measured the saturation magnetiza- 
tion of Ni(<i/vi)/Cu(5 nm) multilayers, for 1 nm< d^i <5 
nm. The data show an extrapolated non-magnetic "dead- 
layer" Ni thickness of 0.25 ± 0.05 nm at each Ni/Cu in- 
terface. Thus we show data only for d^i > 1 nm. In 
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the 1-nm sample one clearly sees the Nb gap, but with 
a significant suppression of the bulk Nb features due to 
the proximity effect in the strong ferromagnet. As we 
increase the Ni thickness, the zero-bias dip in the DOS 
quickly decreases in magnitude. At dNi = 1.5 nm, we 
observe the first sign of an inversion in the differential 
conductance at zero bias, followed by a maximum in- 
version at dNi — 1-75 nm. The features in dl/dV at 
\V\ = A & 1.4 mV have also been inverted but occur 
at the Nb gap voltage in all the samples measured. The 
inversion cycles quickly and by d^i = 2.5 nm, the sam- 
ples return to the non- inverted regime. The dip in the 
DOS reaches its maximum at dNi — 3 nm, then a second 
inversion occurs starting at dNi ~ 4.25 nm. (The second 
inversion is more apparent in the expanded scale of Fig 
2.) This second inverted state looks as though it might 
extend past 5 nm. 



1.025- 
1.020- 
1.015- 
1.010- 
1 .005 ■ 
1.000- 
0.995 ■ 
0.990 ■ 
0.985 ■ 
0.980 ■ 
0.975 ■ 
0.970 ■ 
0.965 ■ 
0.960 ■ 




d Ni ( nm ) 



FIG. 4: (Color online). Black squares represent the normal- 
ized dl/dV at V = for several Ni thicknesses. Red circles 
represent the normalized dl/dV at V — 1.4 mV. The black 
and red lines link the average values of each thickness for 
V = and V = 1.4 mV, respectively. Inset: Data for dNi > 4 
nm on an expanded vertical scale. 

The oscillation is best illustrated in Fig 4 with a plot 
of the normalized differential conductance at both V = 
and V = 1.4 mV versus Ni thickness. The oscillation 
period is irregular; the first inverted region persists for 
dNi ~ 1.3 — 2.4 nm, while the second non- inverted region 
lasts much longer, for dNi ~ 2.4 — 4.3 nm. One can 
see in the inset that the oscillation occurs about a value 
dl/dV ~ 0.9996. This offset in the oscillation is due to 
the aforementioned, weak V-shaped zero-bias anomaly. 

The transitions between normal and inverted DOS 
have been predicted to occur at F-layer thicknesses ex- 
actly half of where 0-tt transitions occur in S/F/S Joseph- 
son junctions. 1 After reviewing the literature on S/Ni/S 
junctions, we see that our junctions are most similar to 
those of Blum et ai 7 and Shelukhin et a/., 9 where a Cu 
buffer layer is placed on each side of the Ni layer. In con- 
trast, Robinson et aZ*ii have Nb in direct contact with 



Ni, with a Ni "dead-layer" thickness of « 0.8 nm at each 
Ni/Nb interface in comparison to our Ni dead-layer thick- 
ness of only 0.25 =b 0.05 nm at each Cu/Ni interface. The 
Ni thickness at which the first 0-tt transition is observed 
by these groups varies quite a bit - 2.6 nm, 1.7 nm (ex- 
trapolated value), and 3.8 nm, for Refs. Q, and (llj . 
respectively. Since our first 0-tt and tt-0 transitions occur 
at 1.3 and 2.5 nm, respectively, we would expect the first 
two transitions in S/Ni/S junctions to occur at 2.6 and 
5 nm, in reasonable agreement with the values observed 
by Blum et ai 7 

Theoretical calculations of the DOS in S/F bilayers 
cover several regimes, defined by the relative strengths 
of the ferromagnetic exchange, E ex , the superconducting 
gap parameter, A, and the impurity scattering, S/r e , as 
well as the relative sizes of dp and the mean free path, 
l e = vpr e . In the dirty limit, the Usadel equations pro- 
vide a clearcut prediction of oscillation of the tunneling 
DOS, with period of order £p = (HD / 'E ex ) x l 2 : , where 
D = vp^e/3. 20 In the clean limit, the predictions are 
less straightforward. Solving the Eilenberger equation in 
the ballistic limit leads to the conclusion that the DOS 
does not oscillate in a semi- infinite ferromagnet j 21 i 22 Os- 
cillations are predicted to occur, however, in the pres- 
ence of weak disorder, with an amplitude proportional 
to h/(E ex r e ) 21 ' 22 With a finite ferromagnet, the oscil- 
lation is regained even in the pure ballistic limit due to 
specular 23 or diffuse 2 ^ scattering from the film bound- 
aries. In constrast, Sun et ai solved the Bogoliubov- 
deGennes equations and claimed that oscillations should 
occur in either a finite or semi- infinite F layer. 25 

Because our observed DOS variations are not peri- 
odic in dNi i we do not attempt to fit our data with Us- 
adel theory. We believe that our samples most closely 
match the assumptions in the papers by Zareyan et 
al. 2A From earlier work on S/Ni systems, we expect that 
E ex « 100 meV)L2J± and v F = 2.8 x 10 5 m/s for Ni^ 
while A = 1.5 meV for Nb. From our measured Ni 
resistivity, pNi = 33 n^m, we deduce l^ 1 = 45 nm. 
This puts our samples in the "intermediate" regime with 
E ex ^> H/r e ^> A. By using Cu buffer layers next to 
the Ni, we limit the scattering events in our junctions 
considerably. As stated earlier, we find a spin diffusion 
length of 21 ±2 nm, low overall spin-scattering asymme- 
try, and very low Cu/Ni interface specific resistance in 
our multilayers: AR Cu / N i = 0.18 ± 0.03/^m 2 . 15 This 
low interface resistance corresponds to a probability of 
scattering of only ~15% at each Ni/Cu interface. 27 The 
Cu/Nb interface is "rough" in the sense that there is 
significant diffusive scattering at this interface, as deter- 
mined from its measured interface specific resistance of 
ARcu/Nb = 1.1 ±0.15/^m 2 . 29 This value of AR Cu /Nb is 
larger than the total AR = 2AR Cu /Ni + PNidNi + PCudcu 
= 0.54 and 0.67 fVtm 2 of the Cu/Ni/Cu region for d m 
= 1 and 5 nm, respectively. We also expect there to be 
diffuse scattering at the CU/AI2O3 (tunneling) interface. 

Plots of the energy dependence of the DOS shown in 
the papers by Zareyan et al£^ agree qualitatively with 
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our data. Performing a quantitative fit of the theory 
to our data, however, is problematic. The theory pre- 
dicts that the first 0-tt transition should occur at very 
small djsn^ a flaw that may be correctable by adding spin- 
dependent interfacial phase shifts to the theory. 30 (One 
could also argue that, because of the 0.25-nm dead layers 
at the two Cu/Ni interfaces, one should subtract 0.5 nm 
from our nominal sample thickness before fitting to the 
Zareyan theory, but that is not nearly enough to bring 
theory into agreement with experiment.) The theory also 
predicts large oscillations in the normalized DOS (i.e., 
large deviations from 1) at zero energy - much larger 
than what we observe in the experiment. The amplitude 
of the theoretical oscillations can be reduced by assuming 
a very small transparency T of the Nb/Cu interface; such 
an assumption, however, is incompatible with the mea- 
sured boundary resistance AR Cu /Nb = 1.1 =t 0.15/^m 2 , 
which implies that T & 0.5. 27 Strong spin-flip scattering 
would also reduce the amplitude of the DOS variations; 
the long measured spin memory length in our Ni films, 
however, precludes that explanation for these samples. 
One could assume that the variation in F-layer thickness 



over the junction area is very large, thereby smearing 
out the oscillations; we believe that such an assumption 
is unrealistic. 

A previous measurement of the Nb/Ni system^ did 
detect signs of one — tt transition, but the data con- 
tained additional low-energy features, which were later 
interpreted as signs of p-wave spin-triplet pairs. 31 We do 
not observe such low-energy features in our data. 

In conclusion, we have observed multiple oscillations 
in the DOS of S/F bilayers as a function of F-layer thick- 
ness, where F is a strong ferromagnet (Ni). The oscil- 
lations can be described qualitatively, but not quantita- 
tively, by the theory of Zareyan et al. 24 Discrepancies 
between theory and experiment may be due to the extra 
Cu layers in our samples, which are not present in the the- 
oretical calculation, or to the absence of spin-dependent 
interfacial phase shifts in the theory. 30 
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